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Stability analysis of the Bianchi type I universe in pure gravity theory is studied in details. We 
first derive the non-redundant field equation of the system by introducing the generalized Bianchi 
type I metric. This non-redundant equation reduces to the Friedmann equation in the isotropic 
limit. It is shown further that any unstable mode of the isotropic perturbation with respect to a 
de Sitter background is also unstable with respect to anisotropic perturbations. Implications to the 
choice of physical theories are discussed in details in this paper. 



PACS numbers: 98.80.Cq; 04.20 -q 



I ■ I. INTRODUCTION 

O ' 

o : 

, Inflationary theory provides an appealing resolution for the the flatness, monopole, and horizon problems of our 
present universe described by the standard big bang cosmology . It is known that our universe is homogeneous and 
isotropic to a very high degree of precision . Such a universe can be described by the well known Friedmann- 
1—9 ■ Robertson- Walker (FRW) metric [§. 

It is also known that gravitational physics should be different from the standard Einstein models near the Planck 
scale For example, quantum gravity or string corrections could lead to interesting cosmological consequences 

Moreover, some investigations have addressed the possibility of deriving inflation from higher order gravitational 
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\q ■ corrections @-[k]]. 

\Q | A general analysis of the stability condition for a variety of pure higher derivative gravity theories is very useful in 

■ many respects. In fact, it was shown that a stability condition should hold for any potential candidate of inflationary 
universe in the flat Friedmann- Robertson- Walker (FRW) space [[Io| . 

In addition, the derivation of the Einstein equations in the presence of higher derivative couplings is known to 
be very complicate. The presence of a scalar field in induced gravity models and dilaton-gravity model makes the 
derivation even more tedious. In order to simplify the complications in the derivation of the field equations, an easier 
way has been described [[TlJJlJ]. We will try to generalize the work in Ref. E3] in order to obtain a general and 
model-independent formula for the non-redundant field equations in the Biachi type I (BI) anisotropic space. This 
equation can be applied to provide an alternative and simplified method to obtain the stability conditions in pure 

■ gravity theories. In fact, this general and model-independent formula for the non-redundant field equations is very 
useful in many area of interests. In particular, it will be applied to study a large class of pure gravity models with 

• *h inflationary BI/FRW solutions in this paper. Any Bianchi type I solution that leads itself to an asymptotic FRW 
KnJ ' metric at time infinity will be referred to as the BI/FRW solution in this paper for convenience. 

Note that there is no particular reason why our universe is initially isotropic to such a high degree of precision. 
Even anisotropy can be smoothed by the proposed inflationary process, it is also interesting to study the stability of 
the FRW space during the post-inflationary epoch. Nonetheless, it is natural to expect that our universe starts out 
as an anisoptropic universe. The universe is then expected to evolve from certain anisotropic universe, e.g. a Bianchi 
type I universe, to an isotropic universe, such as the flat FRW space. Indeed, it was shown that there exists such 
kind of BI/FRW solution for a NS-NS model with a metric field, a dilaton and an axion field p3[ . This inflationary 
solution is also shown to be stable against small field perturbations (m). Note also that stability analysis has been 
studied in various fields of interest Jl5|,[l6| . 

A large class of models with the BI/FRW solutions will be shown to be unstable against arbitrary anisotropic 
perturbations in this paper. We will first derive a stability equation which turns out to be identical to the stability 
equation for the existence of the inflationary de Sitter solution discussed in Ref. Jl^,|l2| . Note that an inflationary de 
Sitter solution in pure gravity models is expected to have one stable mode and one unstable mode for the system to 
undergo inflation with the help of the stable mode. Consequently, the inflationary era will come to an end once the 
unstable mode takes over after a brief period of inflationary expansion. The method developed in Ref. 0,[12| was 
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shown to be a helpful way in choosing physically acceptable model for our universe. Our result indicates, however, 
that the unstable mode will also tamper the stability of the isotropic space. To be more specific, if the model has an 
unstable mode for the de Sitter background perturbation with respect to isotropic perturbation, this unstable mode 
will also be unstable with respect to any anisotropic perturbations. 



II. NON-REDUNDANT FIELD EQUATION AND BIANCHI IDENTITY 

Note that the generalized Bianchi type I (GBI) metric can be read off directly from the following equation: 

ds 2 = g^ehfdx? = -b 2 (t)dt 2 + a 2 {t)dx 2 + aj(t)dy 2 + aj(t)dz 2 , (1) 

with b{t) the lapse function restored on purpose. Note also that the Bianchi type I (BI) metric can be obtained from 
GBI metric by setting the lapse function b(t) equal to one, i.e. b = 1, in equation ([!]). 
Note that one can list all non- vanishing components of the curvature tensor as 

R% = 1 -[H i B + 2B{H i +H*)]5) (2) 
R% l =H i H j B&\ l . (3) 



Here C n kl = ^ m t m ki with e^ k the three space Levi-Civita tensor ||. Here ' denotes differentiation with respect to 
t and Hj = hi/ai is the directional Hubble constant. We have also written B = 1/b 2 for later convenience. Note 
also that the indices i, j in both sides of above equations are open indices. Moreover, R lJ kl = if i = j due to the 
symmetric properties of the curvature tensor. 

Given a pure gravity model one can cast the action of the system as S = J d^x^fgL = 
N J dt{a\a2a^/y/B)L{Hi,Hi,B,B) in the GBI spaces. Here N is a time independent integration constant. We 
will denote the volume factor as V = a\a2a^b. If we take VL as an effective Lagrangian, one can show that the 
variation with respect to b gives 

-~L + ~HiLi + HtL* - i {Hi + 3HH t + H t d/dt)L l = (4) 

after setting B — 1. Here the last term comes from the variation with respect to 5(HiB) with the help of the identity 

' SL -S(H t B) = -J^SiHiB) -> -^-HiSB (5) 
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once we set B = 1. One also needs an integration- by-part with respect to SB in order to obtain the result indicated 
in the last term of Eq. ([l]). Note that the first equality in the Eq. comes from the fact that the factor HiB 
always shows up with 2BHi as indicated by the curvature tensor R tl t j in Eq. (U). Therefore, one ends up with the 
^-equation as 



L - HiU = (Hi - 3HHi - Hid/dt)^ (6) 

after setting B = 1. Here H = Hi/3, Li = SL/8Hi and L 1 = 8L/SHi for convenience. In addition, one can also 
derive, 

L+(3H + d/dt) 2 L i = {3H + d/dt)L h (7) 

as the variational equation of a^. The derivation of this equation is tedious but straightforward. In addition, VL is 
normally referred to as the effective Lagrangian. We will also call L as effective Lagrangian unless confusion occurs. 
Note that the equations (0) are in fact the space-like ij component of the Einstein equation 

G>„ = t^u (8) 

with t^y denoting the generalized energy momentum tensor associated with the system. It is known that one of these 
equations is in fact a redundant equation. Indeed, one can define H^ v = G^„ — and write the field equation as 
H^ = 0. 

Hence one has 



2 



D^H^ = 



(9) 



from the energy conservation (D ll t* lv = 0) and the Bianchi identity (D^G^ V = 0). Now we have three independent 
scale factors Oi and four equations. Therefore, one of the four equations has to be redundant. Indeed, the extended 
Bianchi identity (^|) can be shown to give 

(fit + 3H)H tt + H * H " = °> ( 10 ) 

i 

as soon as the BI metric is substituted into equation (||). Therefore, the equation (|l0| ) indicates that: " H u = implies 
J2i HiHu = 0". Hence, two of the Ha equations vanish will imply that the third one also vanishes. 

On the other hand, H i = H2 = H3 = implies instead (<9 t + 3H)H U — 0. This implies that VH U = constant with 
V = 0,10,20,3. Hence the H t t equation is the non-redundant equation while we are free to ignore one of the three Ha 
equation. Hence any conclusion derived without the H u equation is known to be incomplete. 



III. PERTURBATION AND STABILITY 



One can then apply the perturbation, Hi = HiQ + SHi, to the field equation with Hio the background solution to 
the system. This perturbation will enable one to understand whether the background solution is stable or not. In 
particular, one would like to learn whether a BI/FRW type evolutionary solution is stable or not. It is known that 
our universe could start out anisotropic even evidences indicate that our universe is isotropic to a very high degree of 
precision in the post inflationary era. Therefore, one expects that any physical model should admit a stable BI/FRW 
solution. 

Our result indicates that FRW inflationary solutions with a stable mode and an unstable mode is a negative result 
to our search for a physically acceptable model. Note that FRW inflationary solutions with a stable mode and an 
unstable mode will provide a natural way for the inflationary universe to exist the inflationary phase. Our result 
indicates, however, that such models will also be unstable against the anisotropic perturbations. Therefore, such 
solution will be harmful for the system to settle from BI space to FRW space once the graceful exist process is done. 

First of all, one can show that the first order perturbation equation from the non- redundant field equation (^), with 
Hi — > H + SHi , gives 



L^HSHi 



{3H 2 - H)L 02 + 3HHL 12 + 3HHL 
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+(6HL 01 + 3H Ln - HL 20 + HL 02 + 3HHL 12 + 3HHL 21 )5H = (11) 

where all field variables are understood to be evaluated at the background FRW space where H = H = a/ a for all 
directions, with a the FRW scale factor. We will write H as the Hubble parameter for the FRW space for convenience 
from now on. L a {, = S a+b L/SHi i SHi 2 ■ ■ ■ Hi a SHj 1 SHj 2 ■ ■ ■ SHj b \Hi-*H- I n deriving above equations, we have used the 
following identities 

^ - (13) 
SHi 3SH' [6) 

fC^HiHt)^ f(3H 2 ) (14) 

i 

when we take the limit Hi — > H carrying the system from the BI space to the flat FRW space limit. Here /Q^ H%Hi) 
denotes any functions of the variable Y^i HHi. One can show that the Eq. ([ll]) reduces to 

L 02 SH + 3H L 02 SHi + (6L i + 3H L U - L 20 )SH = (15) 

once we set H = Hq = constant which denotes the de Sitter space. This equation is identical to the stability equation 
for the existence of an inflationary de Sitter solution discussed in Ref. jlOlO]. One notes that there are two other 
independent i/^-equation remain to be checked. These equations can be shown to be redundant after the limiting case 
Hq = constant is implemented. Indeed, the anisotropic perturbation on the Hij -equations are expected to reproduce 
the redundant isotropic perturbation equation shown in Ref. |12] . 



3 



Note that an inflationary de Sitter solution is expected to have one stable mode and one unstable mode for the 
system to undergo inflation with the help of the stable mode. Indeed, one can show that the Eq. ( Jl5| ) can be solved 
to give 

SHi = a exp[B + t] + di exp[-B_i] (16) 
with B± = —3Ho/2 ± y / Ao/2io2 and arbitrary constants Ci,di to be determined by the initial perturbations. Here 



Ao = 9HqLq 2 — 4Lo2(6Loi + 3HoLn — L20) is the discriminant of the characteristic equation for the Eq. (15) 

L 02 x 2 + 3H L a2 x + (6L01 + 3H L n - L w ) = 0. (17) 

There will be an unstable mode if B + > 0. Therefore, the inflationary era will come to an end once the unstable mode 
takes over. It was shown earlier to be a helpful way to select physically acceptable model for our universe. Our result 
shown here indicates, however, that the unstable mode will also tamper the stability of the isotropic space. Indeed, if 
the model has an unstable mode for the de Sitter perturbation, this unstable mode will also be unstable against the 
anisotropic perturbation. 

For example, one can show that the model |l(J 

C = -R-aR% 1 R%R»> llv (18) 

admits an inflationary solution when a < 0. Note that this model is the minimal consistent effective low-energy 
two-loop renormalizable Lagrangian for pure gravity theory fll7[. Indeed, one can show that 



2(H l + H 2 ) 3 -H 



I - 4a(^ H 3 ) 2 -> 6(H + 2H 2 ) - 24a {H + H 2 ) 3 + H 6 (19) 



when we set Hi — > H . Hence one can show that the generalized Friedmann equation (g) gives 

f# = -l/4a. (20) 
In addition, the stability equation ( p"l| ) for SHi can be shown to be 

12aH%5Hi + 36aH$5Hi - (1 + 12aH$)6Hi = (21) 
This equation can be solved to give 



SHi = c,exp[(^35/3- 3)H t/2) + <U exp[-( v / 35/3 + 3)H t/2] (22) 

with arbitrary constants Ci,di to be determined by the initial perturbations. This indicates that this model admits 
one stable mode and one unstable mode following the stability equation ( |l5| ) for the inflationary de Sitter solution. 
It is shown to be a positive sign for an inflationary model that is capable of resolving the graceful exist problem in 
a natural manner. Our result indicates however that this model also admits an unstable mode against anisotropic 
perturbation. Hence this model will have problem remaining isotropic for a long period of time. Therefore, pure 
gravity model of this sort will not solve the graceful exist problem. One will need, for example, the help of certain 
scalar field to end the inflation in a consistent way. 

One expects any unstable mode for a model to be of the form SHi ~ exp[lHot], to the lowest order in Hot, in a 
de Sitter background with I some constant characterizing the stability property of the model. In such models, the 
inflationary phase will only remain stable for a period of the order At ~ 1/IHq. The inflationary phase will start 
to collapse after this period of time. This means that the de Sitter background fails to be a good approximation 
when t 3> Ai. Hence the anisotropy will also grow according to SHi ~ > SHf exp[/_ffoAi] with SH® denoting the initial 
perturbation. 

The measure of anisotropy can be estimated by computing the anisotropy parameter A = yj^j SHf /3H 2 ~ 
x/J2i(3Hi) 2 /3H§ exp[lH t}. Hence A ~ \/J2i(^i) 2 /^Hq exp[60Z] for a typical inflationary model which requires 
a 60 e-fold expansion. This gives us a hope that small anisotropy observed today can be generated by the initial 
inflationary instability for models with appropriate factor /. 

On the other hand, it is known that the Einstein-Hilbert model 

C = -R-2K (23) 

admits only stable modes, which requires SHi — 0, which is bad for the natural graceful exist. This model is, however, 
stable against anisotropic perturbations which tends to keep the universe isotropic as long as the model is in charge. 
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IV. CONCLUSION 



In short, the result of this paper shows that graceful exist and stability of any de Sitter model can not work along 
in a naive way. The physics behind the inflationary de Sitter models appears to be much more complicate than one 
may expect. In another words, the phase transition during and after the inflationary phase deserves more attention 
and requires extraordinary care in order to resolve the problem lying ahead. 

Note that the non-redundant field equation for field theories with many different sorts of fields coupled to the system 
can also be derived similar to the pure gravity models Similar arguments also apply to these theories. Therefore, 
theories with one unstable mode under anisotropic metric perturbations with respect to the de Sitter background will 
not be able to hold the de Sitter space stable for a long period of time even there exist stable modes. Note that, 
in general, one needs to consider two different fourth derivative curvature terms when higher derivative theories are 
considered in four dimensions. 
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